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Abstract 

We study the beam single spin asymmetry in semi-inclusive hadron production in deep inelastic 
scattering in at order 1/Q. There are two competing contributions: the leading order transverse 
momentum dependent parton distribution hj^(x,k±) convoluted with chiral-odd fragmentation 
function e(z), and the chiral-odd distribution function e(x) convoluted with Collins fragmentation 
function H^(z, k'j_). We estimate this asymmetry and compare with the experimental measure- 
ments from CLAS and HERMES collaborations. 
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The single spin asymmetry (SSA) is a novel phenomena in high energy spin physics, 
and has attracted much interest in recent years In particular, the measurements from 
the HERMES, SMC, and JLAB collaborations show a remarkably large SSA in the semi- 
inclusive deep inelastic scattering, such as pion production in 7*0 — > nX, when the proton is 
polarized transversely to the direction of the virtual photon 0,0,0]. On the theoretical side, 
there are many approaches to understanding SSA using Quantum Chromodynamics (QCD) 
phenomenology 0, IE]- Recent interest focuses on the transverse momentum dependent 
(TMD) parton distributions and their implications to the semi-inclusive processes in deep 
inelastic scattering (EBBS EH Hi EH HI- For example, the Sivers function is one of 



those TMD parton distributions representing the asymmetric distribution of quarks in a 
transversely polarized proton, which correlates the quark transverse momentum and the 
proton polarization vector S± [7]. It contributes to the target SSA in semi- inclusive deep 
inelastic scattering. The existence of the Sivers function has been confirmed recently [lj, 

EE III 

Il7l |. where the final state interactions from the gauge link in the gauge invariant 



definition of TMD parton distributions play an important role. 

In this paper, we will study the SSA related to the beam polarization. Unlike the target 
SSA, the beam SSA is subleading in 1/Q, which will eventually vanish as Q 2 — > 00. How- 
ever, at some intermediate Q 2 , but still large enough to guarantee the factorization, this 
asymmetry might be important and measurable. Experimentally, HERMES collaboration 
found this asymmetry consistent with zero 0, but the CLAS collaboration at JLab found 
sizable asymmetry, in the order of 4% E3 • I n the literature 0, 0, 0] , this asymmetry 
has been associated with the twist-3 chiral-odd distribution function e(x) [IE] convoluted 
with the Collins fragmentation function H^-(z,k' ± ) However, this is not the complete 
picture at this order. There is an additional contribution: the leading order TMD parton 
distribution hj;(x, k±) [12] convoluted with a chiral-odd fragmentation function e(z) 
We will demonstrate the existence of this contribution below. 

The TMD parton distribution hj; represents the correlation between the quark's trans- 
verse momentum and polarization in an unpolarized proton state (I2I 22l |. It has the same 



features as the Sivers function: it is a leading-order distribution; it is nonvanishing due 
to the final state interactions; it depends on the quark orbital angular momentum of the 
nucleon j23|. Various model calculations have also shown that it has a similar size as the 
Sivers function (HI, EH IE|- On the other hand, since h± is chiral-odd, it is very difficult 



to probe it in deep inelastic scattering, just like the transversity distribution. Our analysis 
shows that it contributes to the beam SSA, which can be used to extract the distribution 
itself. Another possible way to study hf- distribution is the asymmetry in the Drell-Yan 
process 22l |. 



We first derive the beam SSA in deep inelastic scattering. The semi-inclusive hadron 
production cross section can be expressed as 

** W -L^W , (1) 



dxBdydzd 2 P±h AzxbQ 2 s 

where L^ u and W^ u are leptonic and hadronic tensors, respectively. We work in a frame 
where the virtual photon's momentum q and the proton's momentum P are in the z direction, 
and the incident and outgoing lepton's momenta I and V form a scattering plane. We can 
define the azimuthal angle of any momentum as an angle relative to the scattering plane. The 
variable s is the lepton-hadron total energy square, Q 2 = —q 2 the virtuality of the photon, 
and the dimensionless variables xb, y, and z are defined as xb = Q 2 /2P ■ q, y = P ■ q/ P ■ I, 
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FIG. 1: The relevant diagrams contributing to the beam spin asymmetry associated with h\ <8>e(z) 
for semi-inclusive pion production in deep inelastic scattering. 



z = P ■ p-x/P ■ <i- The variable P±h is the transverse component of the observed pion's 
momentum p n . We further introduce two light-like vectors: p and n, which satisfy p 2 = 0, 
n 2 = 0, p ■ n = 1, p~ = 0, and n + = 0. All momenta can be expressed in terms of p 
and n and the transverse momentum component. The leptonic tensor has the symmetric 
and anti-symmetric parts, = —Q 2 g fiu + 2(Z M ^ + l„l') + 2i\ e e^ vpa l p l' a , where A e is the 
polarization parameter of the lepton. The antisymmetric part will give the beam spin 
asymmetry, convoluted with the antis ymm etric part of the hadronic tensor W£ v . 

We follow the previous studies |9j, LHJ, and use QCD factorization to calculate the 
hadronic tensor W^ v , which can be separated into soft and hard parts. The hard part 
can be calculated in perturbative QCD; and the soft parts are nonperturbative and can 
be parameterized in Lorentz-invariant and gauge-invariant distribution and fragmentation 
functions. The relevant Feynman diagrams are shown in Fig. 1. The lower part of these 
diagrams represent the parton distribution of the target, and the upper part involves the 
fragmentation of a quark into a pion. The symmetric part of W^ v in the leading order has 
a contribution from diagram (a). The antisymmetric part has two contributions. Since the 
contribution from the e(x)®H^-(z, k' ± ) term has been calculated in Ref. [10(, in the following 
we will focus on the contribution from the hj;(x, k±) ® e(z) term. Because e(z) is a twist-3 
fragmentation function, we need to include the diagrams (b) and (c) of Fig. 1 to get an 
electromagnetically gauge-invariant result. 

The parton distributions can be defined from the following density matrix 



M p (x,k ± ) =p H 



di-d\ 



l e -*(?-fc + -f±-fc±) 



(psm~,^cm\ps) 



(2) 



(3) 



where £ is the gauge link The density matrix has the following expansion [l2T ]. 

M p = k ± )j + 2M h U%, k ± )a^k^p v + ... , 

where fi is the usual unpolarized unintegrated parton distribution, and both fi and are 
leading order in twist counting. 

Similarly, for the pion fragmentation functions, we can define the density matrix as 9, 



(4) 



where z = p n ■ p/k' ■ p, and — zk' ± is the transverse momentum of pion relative to the quark's 
momentum. This fragmentation matrix density has the expansion 0|, 



M«=\h{z,k' L )i+\-^- 
2 - 2p w -p 



e(z, k' ± 



+ 



(5) 
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where fi is the usual unpolarized unintegrated fragmentation function; and e is the twist- 
3 chiral-odd fragmentation function. There are two twist-3 and chiral-odd fragmentation 
functions for the pion j2l|, but we only keep the one which contributes to beam SSA. As 
argued in [§, |2l[, instead of the pion mass, we put the nucleon mass as the coefficient in 
front of e(z, k' ± ), because the pion mass vanishes in the chiral limit but the density matrix 
Ai n does not. 

The contribution from Fig. 1(a) to the hadronic tensor W^ u can be calculated as, 

W^ a) = 2z J d 2 k ± d 2 k' ± 5 {2) (P ±h /z -k± + k' ± )ti [M p ^M n Y] ■ (6) 
Substituting the expansions of the density matrices A4 P and A4 n , we have 
w »v{a) = 2z f d 2 k ± d 2 k' ± 5 (2 \P ±h /z-k ± + k' ± )\f 1 (x,k ± )f 1 {z,k' ± ){ P »n u +p»n> 1 -g' 1 ' / ] 



+ih±(x,k ± )e(z,k'J—\p»kl-pW ± ) \ , (7) 

Ptt'P J 

where the electric charge of the quark and the sum over all quark flavor are implicitly 
assumed. The first term in the bracket is the symmetric part of the tensor, and the second 
one antisymmetric. The symmetric part itself is electromagnetic gauge invariant, while the 
antisymmetric part is not, i.e., q^W^ia) ^ o. However, after including the contributions 
from diagrams (b) and (c) in Fig. 1, we can recover the gauge- invariance This leads 
to the following result for the antisymmetric part: i2hj;(x, k±)e(z, k' ± )/z/Q 2 \T^k\_ — T u k±], 
where 7 '' 2xP» • q 1 ' . 

Including also the contributions from the convolution of e(x) with the Collins function 
Hi(z, k' ± ) |lOj], we get the complete result for the antisymmetric part of the hadronic tensor 
at order 1/Q, 

Wr = 2zJ d 2 k ± d 2 k' ± S^ 2 \P ±h /z~k ± + k' ± ) lihi;{x, k ± ) e( - Z ^ A [T»kl - T u k^_] 

-ixe(x, k±)Ht(z, k' x )-^ ]T^l - Tk'f] | , (8) 

and it contributes to the beam SSA. We modified the definition of the Collins function 
Hi in [HI 11 1 by a factor of M p /m w by the same argument we used in Eq. (j3J) for the 
fragmentation function e(z). The definition of e(x) follows |2(|. Since T M is on order of Q, 
we see that the above antisymmetric part will contribute to the cross section in the order of 
1/Q. That means the beam SSA will be 1/Q suppressed, which is different from the Sivers 
effect contribution to the target SSA being the leading order effect. 

Substituting the hadronic tensor into the differential cross section formula Eq. (|I|L we 
will get: 

d 5 a 4ira 2 f - , , , „m,P 



dxBdydzd 2 Pj_h xbV 2 s 



J d 2 k ± d 2 k'^ 2 \^f - k ± + k' ± ) |(1 - y + £)fi{x, k ± )f\(z, k' ± ) 
+ A e — h\ {x, k±)\k±\ sm<p k 

~- K( 2yy// Q — xe ( x > k±)Hi(z, k' ± )\k' ± \ sin0 fc /| , (9) 
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where 4>k and are the azimuthal angles of the momenta k± and k' ± relative to the scattering 
plane, respectively. For example, we define sin^^ = I x I' ■ k±/\l x l'\\k±\. Since only one of 
the two transverse momenta k± and k' ± is relevant for the sin <j) asymmetry, we can integrate 
out the other without assuming the transverse momentum dependence of the distribution 
and fragmentation functions. After that, the <pk or dependence leads to 4>h dependence, 
where <ph is the azimuthal angle of the observed hadron relative to the scattering plane: 
sm<p h = lxl' ■ P± h /\l xl'\\P ±h \. Finally, 

<] ' a f (i-y + y 4) f <?k^\p Lh -zk L )h{xM)h{z) 



dxBdydzd 2 P±h xbV 2 s { 2 

+Ae %V g^ / d2k ^ (2 \ P ^-zk ± )hi(x,k ± )\k ± \^sm<p h (10) 

+Ae %V g^ y^^^^ + ^Dxe^^^fcDI^Isin^.j , 

where a sign has changed in the last term because k' ± and P±h have opposite directions. The 
integrated fragmentation functions are defined as fi(z) = z 2 f d 2 k' ± fi(z, k' ± ), the same for 
e(z), and the distribution function e(x) = f d 2 k±e(x, k±). 

We can further simplify the differential cross section by integrating out the transverse 
momentum P± h but keeping the dependence on (f), 



d 5 a 2a 2 



dxBdydzdcj) xbV 2 s 



{l-y + y 2 /2)f x (x)h{z) 



+ 2\ e y^y^h^ 2 \xf-^^<P 
+2\ e y^T^^xe(x)Ht (1/2 \z) S m^ , (11) 

where the integrated parton distribution fi(x) = J d 2 k_ s _f 1 (x,k_ L ), h^ 1 ' 2 ' (x) = 
J d 2 k ± \k ± \/M p hi(x J k ± ) J and fragmentation H^ (1/2 \z) = z 2 J d 2 k' ± \k' ± \/M p H^( Zj k' ± ). If 



we write the differential cross section as da oc 1 + A y sin i 



T^ l / 2 )(£ 



\ e fd y dzdx B 2 -^^(ht il/2 \x) e -& + xe(x)Ht y ^ 

/ dydzdxB 1 '^^ 2 fi(x)fx(z) 



The xb and z dependence of A y can also be similarly calculated. We note that the two 
contributions have exactly the same dependence on y, which makes it difficult to distinguish 
them experimentally. 

With known distribution and fragmentation functions, we can predict the beam SSA. 
However, up to now, except for the unpolarized quark distribution fi and fragmentation 
fi, these functions can only be estimated in models. In addition, the model calculations 
are not consistent at the present stage. For example, controversial predictions exist for the 
Collins fragmentation function j23, 0|, and we have a wide range of predictions for the 
leading-order TMD parton distribution hj;(x,k±) from models (HI EH EH • So, a reliable 
model prediction is not possible at present. However, we can still gain some insight for these 
functions by comparison with the experimental data. For example, in 0], the beam SSA 



5 





0.1 r 




0.08 


-s- 


0.06 : 


'a) 




ZD 




— I 
< 


0.04 ; 




0.02 



CLAS 



0.5 



0.6 



0.7 



0.8 



FIG. 2: The beam SSA prediction from hj~(x) (g> e(z), compared with the experimental data from 
CLAS 0]. The z dependence solely comes from the ratio of the fragmentation functions e(z) / f\(z). 



has been interpreted as the result of the Collins effect, and the experimental data were used 
to extract the distribution function e(x). 

In this paper, we take an alternative extreme. We interpret the beam SSA as a result 
of the first term in Eq. (|T2j) . To compare with the experimental data, we assume that 
the factorization works at the energy range covered by the experiment. This contribution 
depends on the chiral-odd fragmentation function e(z), which has been calculated in a chiral 
quark model in Ejl . To a good approximation, we have 



e(z 



1 — Z M r 



f 



31 -z 



(13) 



where m q is the constituent quark mass, and M n the nucleon mass. The above relation is 
only true at the scale of A x , and at higher scale their relation might breakdown because 
the evolution of these two functions is different. However, as a rough estimate, we will 
adopt such approximations. The chiral quark model prediction for the usual unpolarized 
fragmentation function f\(z) is consistent with the experimental data after considering the 
evolution effects ji^. We should also note that the chiral quark model is not suitable for 
the calculation of the fragmentation function at z — > 1 region, where the invariant mass of 
the fragmenting quark exceeds the cutoff of the model, A x . 

The z dependence of the asymmetry A y only comes from the ratio of the two fragmen- 
tation functions e(z) and fi(z) in Eq. (|12|) . and the simple relation Eq. (J13)) can be used 
to predict the z dependence of A y . In Fig. 2, we show the normalized asymmetry predic- 
tion from this term compared with CLAS measurements. The most striking observation is 
that this simple relation Eq. ()13|) agrees with the experiment very well. The normalization 
of A y also depends on the TMD parton distribution h^. This distribution involves more 



complicated dynamics 0, |25l 2f|, and hence is less reliable compared to the fragmenta- 
tion function e(z) in Eq. ([13)1 . Nevertheless, from what we have now for h^, we can make 
an order-of-magnitude estimate and compare with experiment. For example, a bag model 
calculation shows that the ratio of h^ l ^ 2 \x) / fi(x) at the kinematic region of the CLAS 
measurement 0.15 < x < 0.4 is about 0.04 for u quark j26[. After taking into account other 
kinematic factors in Eq. (JUJ), the asymmetry A y is predicted to be about 0.05, in rough 
agreement with the CLAS result of 0.038 
crude and the sign is inconsistent. 



although the bag model prediction is very 
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Extending the above estimate to the HERMES kinematics, the beam SSA is at least a 
factor of 5 less than what CLAS has found, consistent with the HERMES measurement • 
This is just the consequence of the beam SSA being 1/Q effect. So, there is no contradiction 
between these two experiments. This also agrees with the observation of 

We note that another interpretation of the beam SSA has been made in (3]]], where the 
photon "Sivers" effect was considered. In Ref. |3(j, an 0{a 2 s ) QCD effect to the beam SSA 
have also been investigated. We did not include these effects in our formalism. 

In conclusion, we have calculated the beam single spin asymmetry in semi-inclusive 
hadron production in deep inelastic scattering. Up to 1/Q, there are two contributions: 
the distribution h±(x,k±) convoluted with fragmentation e(z), and e(x) convoluted with 
Hi(z,k' ± ). A simple chiral quark model prediction of e(z) ~ z/3/(l — z)fi(z) agrees well 
with the experimental data on the z dependence of the asymmetry. Further experimen- 
tal data can provide more information on the extraction of the leading order TMD parton 
distribution 
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